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Low mass neutron stars may be uniquely strong sources of gravitational waves (GW). The neutron 
star crust can support large deformations for low mass stars. This is because of the star's weaker 
gravity. We find maximum ellipticities e (fractional difference in moments of inertia) that are 1000 
times larger, and maximum quadrupole moments Q22 over 100 times larger, for low mass stars than 
for 1.4 Mq neutron stars. Indeed, we calculate that the crust can support an e as large as 0.005 for 
a minimum mass neutron star. A 0.12 Mq star, that is maximally strained and rotating at 100 Hz, 
will produce a characteristic gravitational wave strain of /lo = 2.1 x lO"^** at a distance of 1 kpc. 
The GW detector Advanced LIGO should be sensitive to such objects through out the Milky Way 
Galaxy. A low mass neutron star could be uniquely identified from a large observed spin down rate 
and its discovery would have important implications for General Relativity, supernova mechanisms, 
and possibly nucleosynthesis. 



Rapidly rotating neutron stars can efficiently radi- 
ate gravitational waves (GW) because they involve large 
masses that undergo large accelerations. Laser interfer- 
ometer GW detectors have carried out searches for such 
waves [T]. Although there have been no detections, up- 
per limits have been set. For example, GW radiation 
from the Crab pulsar is constrained to be less than 2% 
of the total spin down energy budget [U [5] . To generate 
GW one needs an asymmetry in the star. This could be 
from a static "mountain" supported by the solid neutron 
star crust. Recently we performed large scale molecu- 
lar dynamics simulations of crust breaking 1^. We find 
that neutron star crust is the strongest material known 
with a breaking stress some 10 billion times larger than 
steel. This is because the long ranged Coulomb interac- 
tions and high pressure suppress many failure modes. As 
a result, the crust can support large mountains, and this 
further motivates GW searches. 

In general, a star with rotational frequency v is ex- 
pected to radiate GW with a frequency 2v. However 
GW radiation at frequency v is also possible from pre- 
cession or from a pinned superfluid component of the 
angular momentum j4j. It can be difficult to search for 
GW from stars with unknown v or in unknown binary or- 
bits. However techniques are being developed to increase 
the sensitivity and reduce the large computational costs 
of these searches |S] . The sensitivity of present detectors 
places important constraints on the quadrupole moment 
(522, distance from earth r, and v of detectable sources. 
For example, for LIGO to have detected the Crab pulsar 
requires an ellipticity, fractional difference in moments of 
inertia e = {I^x — Iyy)/Izz > 10"'' [5]. Ushomirsky et 
al. [B] have calculated the maximum Q22 (and e) that 
the neutron star crust can support. This depends on the 
breaking strain a of the crust which is equal to the break- 
ing stress divided by the shear modulus. Using sa 0.1 
from our MD simulations [3], the Ushomirsky formal- 
ism yields a maximum e ~ 10~^ for a 1.4 Mq neutron 



star, see below. Although 10^^ is larger than previous 
estimates because of the large ct, LIGO is not yet sen- 
sitive to these crustal mountains on the Crab. This is, 
in part, because the Crab spins relatively slowly. LIGO 
is already sensitive to crustal mountains on some more 
rapidly spinning stars. 

The core of a neutron star could involve a high density 
exotic solid such as a spatially nonuniform color super- 
conducting phase [3 IH] • The large shear modulus of this 
solid could support a larger e [SHT^. However this re- 
quires that an exotic solid phase exist at neutron star 
densities and for this phase to become deformed. Pre- 
sumably this would have to happen shortly after the star 
is created in a supernova. 

Alternatively, conventional neutron star crust can sup- 
port larger ellipticities for low mass neutron stars. This 
is because gravity is weaker on low mass stars while the 
strength of the crust is unchanged [T3]. Furthermore, low 
mass neutron stars have very thick crusts where most of 
the star has a nonzero shear modulus and can help sup- 
port mountains. We expect, on very general grounds, 
that low mass self gravitating objects can have large de- 
formations. For example, low mass asteroids can be de- 
formed, while larger mass dwarf planets are nearly spher- 
ical. 

However, it may be difficult to form low mass neu- 
tron stars directly in a supernova explosion, because low 
mass protoneutron stars, that are hot and lepton rich, 
are not gravitationally bound. Perhaps low mass neu- 
tron stars can be formed via fragmentation during pro- 
toneutron star formation or neutron star collisions Il4j . 
If so, low mass stars, formed via fragmentation, could be 
rapidly rotating and deformed. 

We speculate that rapidly rotating low mass neutron 
stars could generate strong gravitational waves. We em- 
phasize that whatever one's theoretical prejudice, the dis- 
tribution of neutron star masses is an important obser- 
vational question. Gravitational wave observations could 
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provide unique information on low mass stars. The struc- 
ture of low mass stars is interesting because their cen- 
tral densities are comparable to nuclear density. As a 
result, their structure is closely related to properties of 
conventional atomic nuclei [TS]. Low mass stars should 
have maximum rotation frequencies lower than the over 
1 kHz frequency expected for a 1.4 Mq star. Low mass 
stars, spinning close to their maximum rate, may radiate 
gravitational waves with frequencies near 200 Hz. This 
is where present ground based interferometers are most 
sensitive. 

We describe the structure of low mass neutron stars 
using the equation of state (EOS) of Douchin and Haensel 
for the inner crust and core [16 and the original Baym, 
Pethick and Sutherland (BPS) EOS for the outer crust 
[T7] . We discuss the sensitivity of our results to the EOS 
below. The angle averaged shear modulus ^ of the crust 
is taken from our recent molecular dynamics (MD) results 

m, 

^ = 0.1106 n,. (1) 

a 

Here the ions have charge Z, density rii and the ion sphere 
radius is a = 

(3/47rn,)i/3^ This result for ^ is about 10% 
smaller than the earlier Monte Carlo results of Ogata et 
al. because they neglected electron screening. 



the crust are plotted. Note that Z and A change contin- 
uously in the inner crust because of the statistical model 
empolyed for the EOS [16]. Also shown in Fig. [l]is the 
ratio of the shear modulus to the baryon density ^i/ Pb- 
The star has a radius of 29 km and is mostly solid crust 
with only a small liquid core. For comparison, the struc- 
ture of a conventional 1.4 Mq star has only a thin solid 
crust and a large liquid core. Figure [2] plots the radius of 
a neutron star versus its mass. Low mass neutron stars 
can be much larger than 1.4 Mq stars. 
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FIG. 1: (Color on line) Profile of a low mass 0.12 Mq neu- 
tron star. The charge Z (red dashed line) and mass number 
A (green dashed-dotted line) of nuclei in the crust are plot- 
ted versus radius r. Note that A is arbitrarily set to one in 
the core. The thick line shows the baryon density in fm~^ 
multiplied by 2000 for clarity. The dotted blue line shows the 
mass fraction of free neutrons Xn in the inner crust. Finally 
the dashed two dots line shows the ratio of shear modulus to 
baryon density ^J./ pb in MeV. 

The structure of a 0.12 Mq neutron star is shown in 
Fig. [l] The charge Z and mass number A of nuclei in 



FIG. 2: (Color on line) Radius r (solid black curve and left 
hand scale) and moment of inertia I (dashed red curve and 
right hand scale) versus mass M of neutron stars. 

The mass quadrupole moment Q22 of the star is, 

Q22 = 2j dVp(r)r2Rer22, (2) 

where p(r) is the (energy) density and 122 is a spheri- 
cal harmonic [S]. Ushomirsky et al. have estimated the 
maximum Q22 that can be supported by a crust with 
stress tensor In a related calculation, Haensel et 

al. calculated how the energy of a star rises with de- 
formation |20| . Knippel and Sedrakian assume the com- 
ponents of tij are stressed just to the breaking point in 
a way that maximizes Q22 [H], irr — 2(327r/15)^/^/i(T, 
tr± = 2(167r/5)i/VCT, tA = 6(3271/15)1/^^. Note that 
these relations differ by a factor of two from those in ref. 
[6]. We assume the breaking strain a is, 

a = 0.1, (3) 

based on our large scale MD simulations This gives 
a maximum Q22 of [5], 

(4) 
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Here g{r) is the acceleration due to gravity and U = 
2 + dlng/d\nr. The eUipticity e is the fractional difference 
in moments of inertia, 



d away, is 



1/2 Q22 



(5) 



where we have used Eq. [2] for Q22- The moment of 
inertia / = Izz, calculated in a slow rotation approxima- 
tion 1^T\, is shown in Fig. [2j Note that we neglect the 
effects of rapid rotation which can significantly distort 
the structure of low mass stars because of their extended 
crust ^13j . The moment of inertia, of very low mass stars 
increases because of their large radius. 

The maximum Q22 and e that the crust can support 
are plotted in Fig. |3] The maximum e for low mass stars 
can be as large as 5 x 10^"^. The crust can support a Q22 
that is over 100 times larger, and an e that is 1000 times 
larger, for a low mass star than for a I.4 Mq star. 
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FIG. 3: (Color on line) Maximum eUipticity e = e (solid curve 
and right hand scale) and Q22 (dashed red curve and left hand 
hand scale) versus mass M. 

The maximum rotation frequency J^max, for a star of 
mass M and radius r, has been estimated to be [22l 
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(6) 



If this equation can be extrapolated to very low masses, 
it gives t'max ~ 70 Hz for the O.I2M0, r = 29.5 km star 
in Fig. [1] Explicit two dimensional calculations of the 
structure of a rapidly rotating 0.13 Mq star by Haensel 
et al. give i/^ax = 100 Hz using the SLy EOS [B]. This 
is only slightly larger than the above estimate. Note that 
this star rotating at 100 Hz is significantly oblate with 
the ratio of polar to equatorial radii Tpoic/fcq = 0.7 [T^ . 

A characteristic gravitational wave strain Hq, from a 
star with e and / spinning at frequency v and a distiance 



(7) 



For a O.I2M0 maximally stressed star, spinning at 100 
Hz and at a distance of 1 kpc, this gives a relatively 
strong signal. 



2.1 X 10 
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The LIGO all sky search [24] may have already ruled 
out such an isolated star within distances of order 1 to 
3 kpc, depending on orientation. However the sources 
rapid spin down rate or presence in a binary system could 
be complications that might greatly reduce the excluded 
distance. Advanced LIGO should be sensitive to these 
stars throughout the Milky Way Galaxy. 

Strong gravitational wave radiation will lead to a rapid 
spin down. Therefore, one may be most interested in rel- 
atively young objects. For example, Imshennik presents 
a speculative scenario for SN1987a that involves a very 
rapidly rotating core that fragments to produce a low 
mass neutron star [SFj. We note that GW from a low 
mass star could be strong enough to be detected, even 
given the relatively large distance to SN1987a. There- 
fore GW observations of SN1987a can probe the possible 
formation of a low mass neutron star. 

We discuss the sensitivity of our results to the EOS and 
approximations that we use. First, we assume cold cat- 
alyzed matter. In principle, the composition of accreted 
matter could be significantly different. Haskell et al. have 
studied mountains on normal mass stars made of both ac- 
creted and non-accreted matter |26) . However, it seems 
difficult to make low mass neutron stars out of accreted 
matter. Second, there could be sensitivity to the EOS of 
cold catalyzed matter. Low mass neutron stars have rela- 
tively low densities, and there may be only a small range 
in pressure for realistic EOS at these densities. However, 
the model dependence of the shear modulus, Eq. [T] and 
the breaking strain, may be more important. The crust 
is strongest (and can make the largest contribution to 
the maximum Q22) at high densities near the base of the 
inner crust. Therefore, one is sensitive to the composi- 
tion Z, A of nuclei deep in the inner crust. This com- 
position is set by the symmetry energy and there can be 
significant variations among different EOS models. For 
example, Steiner and Watts find large variations in the 
speed of shear waves deep in the inner crust [27] . In addi- 
tion, we expect non-spherical nuclear pasta phases at the 
base of the inner crust, see for example [5S]. The impact 
of these phases on both the shear modulus and breaking 
strain should be studied further in future work. Finally, 
we use the formalism of Ushomirsky et al. to estimate the 
maximum Q22 ■ This involves the Cowling approximation 
and may omit some boundary terms !26]. Furthermore 
we assume a nearly spherical star, while rapidly rotating 
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low mass stars can be significantly oblate. These approx- 
imations should be improved in future work. Finally, and 
perhaps most difficult, one should study formation mech- 
anisms for low mass stars and ways that the star could 
become deformed. 

We now speculate on the strongest possible sources 
of continuous gravitational waves for ground based de- 
tectors. Objects significantly larger than neutron stars 
can not spin fast enough to radiate at frequencies where 
ground based detectors are most sensitive, see Eq. |6] 
Sources involving static deformations of conventional 
matter are constrained by the breaking stress t « /zct of 
neutron star crust, see Eqs. |1|3[ Note that neutron star 
crust is the strongest material known. Given the source's 
size and composition, we now optimize its mass. Figure 
[3] clearly shows that the smaller the mass, the stronger 
the possible gravitational wave source. We expect that 
this is a very general result. Gravity is weaker for smaller 
masses. This allows the same material strength to sup- 
port larger deformations. However there is a lower limit 
to neutron star masses, if they are to be bound by grav- 
ity. The minimum neutron star mass for the EOS we use 
is 0.094 Mq. We conclude that low mass neutron stars 
could be uniquely strong gravitational wave sources. 

There may be a useful analogy between gravitational 
wave searches for low mass neutron stars, and doppler 
shift searches for hot jupiter exoplanets. Originally some 
astronomers thought that hot jupiters would be unlikely 
to form. However, many planet searches are strongly bi- 
ased towards the discovery of hot jupiters, because they 
produce large doppler shifts. It is important to consider 
this bias in interpreting planet search results. Likewise, 
gravitational wave searches are sensitive to deformation. 
Therefore, they may be strongly biased towards the dis- 
covery of low mass neutron stars, because these stars 
can have large deformations. One should consider this 
bias when interpreting the results of gravitational wave 
searches. 

To actually have a low mass gravitational wave source 
requires three things. First, low mass neutron stars must 
be produced. This is problematic because the minimum 
mass of a hot protoneutron star is much larger. Perhaps 
low mass stars can be formed by fragmentation. Second, 
the star must be deformed. This deformation could have 
occurred when the star was formed in a violent event. Fi- 
nally, and perhaps more straight forward, the star must 
be rapidly rotating. Again if the star was formed via frag- 
mentation it is possible that the star was formed rapidly 
rotating. However, the star will quickly spin down be- 
cause of gravitational wave radiation. Therefore, the sys- 
tem may have to be young in order for it to be still rapidly 
rotating. 

We strongly encourage searches for gravitational waves 
from low mass neutron stars. These stars could be very 
strong GW sources and can be uniquely identified be- 
cause of their rapid spin down rate, that is inversely pro- 



portional to the moment of inertia /. Searches should 
be extended to include large spin down rates. Figure 
[2] shows that / can be an order of magnitude or more 
smaller than that for a 1.4 Mq star. A small observed 
/ would provide a "smoking gun" that the GW source 
is a low mass neutron star. This star could possibly be 
made of exotic matter. However the identification of low 
mass is unique because a larger mass object would likely 
collapse to a black hole before its / became this small. 

The GW detection of a low mass NS would have dra- 
matic implications for General Relativity, Supernovae, 
and perhaps nucleosynthesis. By proving that low mass 
stars are formed, presumably be fragmentation, this dis- 
covery would stress the role of rotation in at least some 
supernovae. This is important because our understand- 
ing of supernova mechanisms may be significantly in- 
complete. Furthermore by demonstrating the ejection of 
large amounts of neutron rich matter, the formation of 
low mass neutron stars may be related to rapid neutron 
capture nucleosynthesis of heavy elements. 
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